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We test the AdS/CFT correspondence by calculating Wilson loops in JV = 4 super Yang-Mills 
theory on R x S 3 in the planar limit. Our method is based on a novel large-rV reduction, which re- 
duces the problem to Monte Carlo calculations in the plane-wave matrix model or the BMN matrix 
model, which is a Id gauge theory with 16 supercharges. By using the gauge-fixed momentum- 
space simulation, we obtain results respecting 16 supersymmetries. We report on the Monte Carlo 
results for the BPS circular Wilson loop, which reproduce the exact result up to strong coupling. 
As a future prospect, we calculate a track-shaped Wilson loop from the gravity side, which shows 
that a clear test of the AdS/CFT for the non-BPS case is also feasible. 
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1. Introduction 

The gauge-gravity duality has been one of the most important subjects in string theory over 
the past decade. The most typical example is the so-called AdS/CFT correspondence between type 
IIB superstring theory on AdS 5 x S 5 and 4d JV = 4 SU(A0 super Yang-Mills theory (SYM). Even 
in this case, however, a complete proof of the duality is still missing. One of the reasons is that the 
parameter region described by the classical supergravity on the string theory side corresponds to 
the strongly coupled region in the planar large-A^ limit on the gauge theory side. In order to study 
the strongly coupled 4d jV = 4 SYM from first principles, one needs to have a non-perturbative 
formulation such as the lattice QCD. The problem here is that the lattice regularization necessarily 
breaks translational symmetry, which is included in the supersymmetry (SUSY). In order to restore 
SUSY in the continuum limit, one generally has to fine-tune parameters in the lattice action 1 . 

In this work we adopt an alternative regularization method based on the idea of the large- 
N reduction [Q], which preserves 16 supersymmetries. Since 4d jV = 4 SYM has conformal 
symmetry, the theory on R 4 is equivalent to the theory on R x S 3 through conformal mapping. 
The novel large-Af reduction relates this theory to a reduced model, which can be obtained by 
shrinking the S 3 to a point. The resulting one-dimensional gauge theory with 16 supercharges can 
be studied by using the gauge-fixed momentum-space simulation [||] as in recent studies of the D0- 
brane system [g]. Thus we can perform Monte Carlo calculations in 4d JV = 4 SYM respecting 
SUSY maximally and without fine-tuning 2 . 

We are going to test the AdS/CFT correspondence by calculating BPS and non-BPS Wilson 
loops in 4d JV = 4 SYM 3 . In particular, we reproduce an exact result for the circular Wilson loop, 
which serves as a check of our method. For a non-trivial test of the AdS/CFT correspondence, 
we consider the track-shaped Wilson loop as an example of non-BPS operators, which cannot be 
calculated by analytic methods relying on SUSY. We calculate it on the gravity side by numerically 
solving a classical string equation of motion. 



2. Large-N reduction for jV = 4 SYM on R x S 3 



Let us first discuss the novel large-Af reduction for Jf = 4 SYM on R x S 3 . By collapsing the 
S 3 to a point, we obtain the plane wave matrix model (PWMM) or the BMN matrix model [|To|] 4 , 
whose action is given by 



Spw 



<?pw 



dt tr 



{D t X M f - - [X M ,X N ] 2 + ~V*D t V - -^Ym[X m ^} 



+^{Xif + y(^) 2 + imjkXiXjXu + i^L^y m y 



(2.1) 



?3 „„ d —2 
I 1 ' 

by D t = d t — i[A, • ], where A(t), as well as X^it) and *P(f)> is an ,/V x N hermitian matrix. The 



Here the parameter \x is related to the radius of S as R S 3 = ~, and the covariant derivative is defined 



'Recently it was claimed that fine-tuning can be avoided at least up to the 1-loop level by using a lattice formulation 
with topological twist 

2 See refs. [7] for proposals for finite N. 

3 See refs. [8, ^] for some preliminary results on the Wilson loop and correlation functions. 

4 Properties of this model at finite temperature are studied at weak coupling JTT| , |T^ ] and at strong coupling Jl3||. 
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range of indices is given by 1 < M,N < 9, 1 < i, j, k < 3 and 4 < a < 9. The model has the SU(2|4) 
symmetry with 16 supercharges. 

The PWMM possesses many discrete vacua representing multi fuzzy spheres, which are given 
explicitly by 

V v 

Xi = M0(l.^®lfe) with £n I k I =N, (2.2) 

7=1 /=1 

where L ; - are the r-dimensional irreducible representation of the SU(2) algebra [Lj r \l!p] = 

(r) 

iSjjkL k . These vacua preserve the SU(2|4) symmetry, and are all degenerate. 

In order to retrieve the planar jV = 4 SYM on R x S 3 , one has to pick up a particular back- 
ground from (2.2), and consider the theory (2.1) around it. Let us consider the vacuum defined 
by 

v + 1 

k I = k, m = n+I — for /= l,--- ,v , (2.3) 

and take the large-/V limit in such a way that 

Tl J?^ k 

/c-^oo, >«, y-^oo, with Apw = PW fixed . (2.4) 

V n 

Then the resulting theory is claimed [Q] to be equivalent 5 to the planar limit of ,jV = 4 SYM on 
Rx S 3 with the 't Hooft coupling constant given by 

^SYM = 27T 2 /lpwC%0 3 = ' (2.5) 

The above equivalence may be viewed as an extension of the large-/V reduction [^J, which as- 
serts that the large-A^ gauge theories can be studied by dimensionally reduced models. The original 
idea for theories compactified on a torus can fail due to the instability of the U(1) D symmetric vac- 
uum of the reduced model JT^]. This problem is avoided in the novel proposal since the PWMM 
is a massive theory and the vacuum preserves the maximal SUSY. This regularization respects 16 
supersymmetries, which is half of the full superconformal symmetry of jV = 4 SYM on Rx S 3 . 
Since any kind of UV regularization breaks the conformal symmetry, this regularization is optimal 
from the viewpoint of preserving SUSY. 



3. Wilson loops in 4d jV = 4 SYM 

Let us consider the following type of Wilson loop 

W (C) = -^Tr & exp j> ds (iA^x^ (s) + \x» (s) | Xjf * fl ) , (3.1) 

where x^(s) = dx ^^ and d a is a constant which satisfies 6 a 6 a = \. The fields A^f andXf 4 represent 
the gauge field and the six scalars in 4d jY = 4 SYM on R 4 , respectively. Due to the particular way 

5 See refs. jl4] | for earlier studies that led to this proposal. This equivalence was checked at finite temperature in the 
weak coupling regime [113] . It has also been extended to general group manifolds and coset spaces |Ea|. 
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in which the scalars appear, one can obtain predictions from the gravity side based on the AdS/CFT 
correspondence as 



lim 

°,/l-SYM- 



W(C) 



-m 



(3.2) 



iV->°°,AsYM->°° \ ' ' SYM 

where S(C) represents the area of the minimal surface spanning the loop C on the boundary of the 
AdS space. 

For the circular Wilson loop W(C c i rc ), which is a (l/2-)BPS operator, there is an exact result 
on the gauge theory side, which is obtained by summing up planar ladder diagrams or by using the 



localization method [17]. The result is given by 



lim 

N— ¥<" 



W(C C - 



SYM 



2k 



■SYM 



f) 1/2 (2A SYM ) 3 / 4 



for k> 



SYM 



»1, 



(3.3) 



(3.4) 



where I\ (x) is the modified Bessel function of the first kind. The result is independent on the radius 
of the circle, which is a consequence of the scale invariance of jV = 4 SYM. At strong coupling 
it agrees with the result obtained from the dual geometry >S(C c i rc ) = — \/2Asym [18]. This is an 



explicit example of the AdS/CFT correspondence. We use the exact result (3_3) for arbitrary Asym 
to check our calculation method. 




1.8 



const.+anti-parallel lines 

circular O 




3.5 



Figure 1: (Left) The track-shaped Wilson loop has two straight lines of length T, which are anti -parallel to 
each other, and two semi-circles with radius L/2. (Right) The numerical results for ' Klck(rx> is plotted as 

V^SYM 

a function of T/L. The circle at T — represents the circular limit (T —> 0) and the straight line represents 
a fit to the anti-parallel-line limit (T — > °°) with an additive constant treated as a free parameter. 

The track-shaped Wilson loop has two parameters T and L as depicted in figure |l| (Left). It 
has two limits: the T — > limit, which corresponds to the circular Wilson loop, and T — > °°, which 
corresponds to two anti-parallel lines. Since the track-shaped Wilson loop for T ^ is a non- 
BPS operator, it cannot be calculated on the gauge theory side by such analytic methods as the 
localization method which rely on SUSY. The limit (T — > oo) of two anti-parallel lines is calculated 



based on the AdS/CFT correspondence from the gravity side as [19] 

1 . / , , \ nr. 47T 2 



lim -ln(W(C track(rx) } 

T \ I SYM 



■V2k: 



1 



■SYM 



r(i/4) 4 L 



(3.5) 



4 



Testing AdS/CFT by Monte Carlo calculation of Wilson loops 



Masazumi Honda 



which is consistent with the conformal symmetry of the theory, and also with the fact that jV = 4 
S YM on R 4 is in the Coulomb phase. Here we extend the calculation on the gravity side to arbitrary 
T/L, and obtain explicit numerical results by solving the classical string equation of motion. In 
figure [I] (Right) we plot the minimal surface obtained by Newton's method. 



4. Monte Carlo method 



In order to simulate the PWMM (2.1), we compactify the ^-direction to a circle of circumfer- 
ence j3. Since we are interested in the properties at zero temperature, we impose periodic boundary 
conditions on both scalars X/\f(t) and fermions *Pa(0> which keep SUSY intact. In Fourier-mode 
simulation [g], we first fix the gauge symmetry completely by choosing A(t) = ^diag(ai , • • • , a#) 
with —% < a a < 71, and then make a Fourier expansion Xj^it) = Y 1 n=-A^M,n^" ont = jf) and 
similarly for the fermions. The upper bound A on the Fourier modes plays the role of the UV 
cutoff. The original PWMM can be retrieved by just taking the limits j3 — > °° and 4 — > oo since 
there are neither UV nor IR divergences. The model regularized by finite j3 and A can be simulated 
by the RHMC algorithm. This method has been applied extensively to the DO-brane system cor- 
responding to n = 0, and the results confirmed the gauge/gravity duality for various observables 
[^]. 6 Since the parameter jx in the action (JO]) can be scaled away by appropriate redefinition of 
fields and parameters, we take \i=2 (R S 3 = 1) without loss of generality. 

The Wilson loop in JV = 4 SYM can be calculated in PWMM in the following way. When we 
perform the conformal mapping from R 4 to R x S 3 , the radial and angular directions are mapped to 
the time and 5 3 -directions, respectively. Therefore, an arbitrary loop on a plane in R 4 is mapped 
to a loop on R x S 3 , which can be projected to a great circle on S 3 . Such a Wilson loop can be 
represented in the large-Af reduced model as 

W red (C) = ^Tr 9 exp £ ds (iA £ + iX^ (s) + |jfc„ (s) \X a d)j , (4.1) 

where e'j = e'j(x^(s)) is the dreibein on S 3 . The expectation value of this operator is related to the 
average of the original Wilson loop as [ 21 ] 



W(C)) =(W red (C)), (4.2) 

/SYM \ / 

where (•••) on the right-hand side denotes the expectation value in the large-Af reduced model 
(PWMM). In the case of circular Wilson loop, the relation ([O]) was confirmed by reproducing 



the SYM result (3.3) from the reduced model to all orders in perturbation theory assuming that 



non-ladder diagrams do not contribute [p 



5. Numerical results 

In figure ^| we present our preliminary results for the circular Wilson loop. We have performed 
the A — > oo extrapolation using A = 6,8,10,12 assuming that finite A effects are 0(1/A). The 
background is chosen to be (n, v) = (|,2) and we performed an extrapolation to k = oo using the 

6 See refs. [E(l for Monte Carlo calculations based on the lattice regularization. 
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data for k = 2,3,4,5 assuming that the finite-k effects are 0(\/k 2 ). We also plot the exact result 
(3.3). Except for the data point at \/Asym = 4, the agreement with the exact result is promising. 
Note, in particular, that we already start to observe a bent from the weak coupling behavior towards 
the strong coupling behavior. This is remarkable considering the rather small matrix size. We 
consider this as a result of the fact that our formulation respects sixteen supersymmetries. 
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Figure 2: The log of the circular Wilson loop normalized by \Msym is plotted against \Msym- The values of 
gpy/N and j3 (and the corresponding values of Asym) we use are listed in the table. The solid line represents 
the exact result (3.3). The dashed line represents the behavior (3.4) at strong coupling, whereas the dotted 
line represents the leading perturbative behavior In (W) ~ zAsym- 



6. Summary and discussions 

We have investigated nonperturbative properties of the 4d SU(°°) JV = 4 SYM from first- 
principles respecting 16 supersymmetries. In particular, we have reproduced the exact result ( [Q| ) 
up to Asym — 10.4. As a nontrivial check of the AdS/CFT correspondence, we are planning to study 
the track-shaped Wilson loop. We have obtained explicit results on the gravity side by numerically 
solving a classical string equation of motion. The results nicely interpolate the two limits, which 
are already calculated, i.e., the circular loop and the anti -parallel lines. We hope to report on the 
calculations on the gauge theory side in the forth-coming publication. 
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